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DESCRIPTION AND ANALYSIS OF A ROCKET-VEHICLE 
EXPERIMENT ON FLUTTER, INVOLVING WING 
DEFORMATION AND BODY MOTIONS 
By H. J. Cunningham and R. R. Lunds trom 

SUMMARY 


Flight tests and a mathematical analysis were carried out to demon- 
strate and confirm a type of subsonic flutter involving rigid-body 
motions and wing deformations. For the configuration considered the 
period of the oscillation was approximately 100 chords per cycle which 
is well within the range of period found in dynamic-stability work on 
rigid aircraft with free controls. A mathematical analysis based on 
two-dimensional incompressible flow provided a conservative prediction 
of the airspeed at which the low-frequency flutter occurred. It was 
found that wing bending stiffness is the important parameter for pre- 
venting such flutter. 


INTRODUCTION 


Interaction of deformations of an aircraft structure with the- passing 
air stream can lead to the dynamic instability known as flutter. For 
bending-torsion wing flutter the frequency of oscillation is fairly high, 
usually approaching the natural torsional frequency of the wing in still 
air. Such an oscillation may be contrasted with ordinary dynamic- 
stability phenomena involving rigid-body modes which lead to much lower 
frequencies that are usually controllable. 

The fact that the calculated flutter speed may be modified by the 
addition of free-body modes has been recognized for many years. For 
example, about 20 years ago it was found analytically in reference 1 
that body mobility had a slight favorable effect on the calculated 
flutter speed of a particular configuration typical of that day. The 
problem has from time to time been reconsidered in both American and 
British literature and the necessity for determining any potential 
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detrimental effect of special configurations, sweptback wings, higher 
speeds, and higher altitudes has lately become more insistent. A recent 
paper by Broadbent (reference 2) considers chiefly by calculation (with 
the mention of some preliminary qualitative experimental results) the 
necessity for including free-body modes in the study of sweptback wings. 

Controlled experimentation involving free-body modes is highly desir- 
able though difficult even at low speeds in a wind tunnel. In reference -3 
Lambourne gives some experimental results and' describes various diffi- 
culties encountered, principally the difficulty of supporting the model 
so that actual flight behavior is sufficiently well simulated. The 
National Advisory Committee for Aeronautics has been engaged for several 
years in making flatter experiments in the transonic range by the use of 
modern telemetering techniques with rocket-powered vehicles. Some of the 
results of these experiments are given in several NACA .papers including 
references k to 6. Such flutter research has been concerned almost 
entirely with wing bending-torsion flutter. Recently, however, in con- 
nection with some, of these experiments flutter failures were obtained 
which definitely involved much lower frequencies than those obtained in 
ordinary bending-torsion flutter. A preliminary account of an unexpected 
low-frequency failure is given in reference U, and two other low-frequency 
failures are reported *in reference 5. 

It was decided to repeat the experiment of the D model of reference U 
with more extensive instrumentation in order to obtain information specific 
to this type of flutter. The present paper includes a description of and 
results from model D of reference U as well as similar and more compre- 
hensive material for the repeat experiment which has been designated 
model E. Results of an analysis developed on the basis of certain sim- 
plifying assumptions are given, and a comparison is made between analyt- 
ical and experimental results. 


SYMBOLS 


a + Xa 


A 


g 


b 


nondimensional position of wing elastic axis ( 

/ \ 100 

nondimensional position of wing-section center of gravity 



geometric aspect ratio of one exposed wing panel, for 
rectangular plan-form wing's (Length/Chord) 

wing semichord in feet 
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El 

s 


^h2 

ft 

fa 


Eh 

Sa 

g 

GJ 

k 

l 

m 

M 

M cr 


bending rigidity, pound- inches^ 

nondimensional distance (in wing semichords) from wing 
midchord line to missile center of gravity measured 
parallel to the missile longitudinal axis, positive 
when center of gravity is behind wing midchord line 

first bending natural frequency, cycles per second 


second bending natural frequency, cycles per second 
first torsion natural frequency, cycles per second 
uncoupled first torsion frequency relative to elastic axis. 


cycles per second ft 


1 - 


( x a/ r q) ' 


- (V0‘ 


1/2N 


structural damping coefficient in bending 

structural damping coefficient in torsion 

mathematical quantity having the qualities of a structural 
damping coefficient 

torsional rigidity, pound- inches^ 


reduced frequency parameter (cob/v); n/k is the period of 
the oscillation in chords per cycle 

length of wing along leading edge outboard of body, inches 

mass of wing per unit length, slugs per foot 

Mach number 

theoretical Mach number at which sonic velocity is first 
attained over section of wing taken perpendicular to 
leading edge at zero lift 


nondimensional radius of gyration of wing section about 

elastic axis (^I a /mb 2 ) where I a is the mass moment 

of inertia of the wing about its elastic axis per unit 
length in slug-feet 2 per foot 
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v 

bco h 

x cg 

x ea 


k (stnd.) 
P 


CO 


COh 


o>a 


a 


time after rocket launching, seconds 
speed, feet per second 

nondimensional flutter-speed coefficient 

distance of center of gravity of wing section behind 
leading edge, percent chord 

distance of wing elastic axis behind leading edge, 
percent chord 

mass ratio (npb^/m) 

value of k when p is standard sea-level density 

air density, -slugs per cubic foot 

angular frequency of vibration, radians per second 

angular uncoupled first bending frequency, radians 
per second 

angular uncoupled first torsion frequency about elastic 
axis, radians per second (co a = 2nf a ) 


■flutter parameter 


* *>) 


EXPERIMENTAL INVESTIGATION 
Apparatus and Methods 


The test vehicle was essentially a. tailless configuration having 
the test wings as the only stabilizing surfaces in pitch. The model was 
powered by a modified Aerojet 12AS-1000 D rocket motor capable of 
carrying it to greater- than-sonic speed with an acceleration of about 
U times gravity. A sketch of model E is shown in figure 1. Center of 
gravity, weight, and moment of inertia in pitch varied with flight time 
because of the fuel consumed. These parameters plotted against flight 
time for model E are shown in figure 2. The launching was made from a 
near-zero-length launcher, as described in reference 6. - Photographs of 
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model E in launching position are shown in figure 3. The test wings 
were made of laminated white pine with an inlay of 0.032-inch 
2I4ST aluminum alloy to duplicate as closely as possible the test wings 
of model D of reference U. Strain gages were mounted on the elastic 
axes of the wings U inches from the wing root, as shown in figure 1. 

A six-channel telemeter was installed in the model with instrumen- 
tation to give continuous readings of the strains, left-wing bending, 
right-wing bending, and right-wing torsion, at the corresponding gage 
locations, and also to give continuous readings of longitudinal and 
normal accelerations of the model center of gravity and angle of attack 
of the missile. Speed of the model was obtained by integration of longi- 
tudinal acceleration, and the altitude was obtained from a pulse-type 
tracking radar unit. Atmospheric conditions prevailing at the time of 
flight were obtained from a radiosonde. 

The accuracy of the quantities measured is believed to be within 
the following limits: longitudinal and normal accelerations, 2 percent.; 

velocity, I4 percent; wing bending and twisting moments, 20 percent; and 
phase angle between any two quantities, 15°. 


Experimental Results 

The physical characteristics of the wings as determined by preflight 
ground tests are listed in table I. Data for model E are presented in 
figures It to 7 as functions of the flight time. Included in these fig- 
ures are the available data for model D from reference U for comparison. 
Figure Ij. includes flight velocities; figure 5 shows Mach number and air 
density; figure 6 shows longitudinal acceleration. Figure 7 presents 
the following data: bending moments of both wings (positive, wing bent 

down) at the strain-gage locations, twisting moment of the right wing 
about its elastic axis (positive, leading edge up) at the strain-gage 
location, and normal acceleration of the missile center of gravity 
(positive, up). It was desired to obtain the angle of attack of model E 
directly as a function of time, but the angle-of -attack indicator was 
inoperative during flight. As a result, the angle of attack, was deter-, 
mined by use of the 'normal acceleration and an assumed lift-curve slope 
of 0.08 per degree as obtained from lifting-surface theory. 

Comparison of figure 7 for model E and figure 9 of reference 1* for 
model D indicates that both flights had quite similar behavior. At a 
speed of about UOO feet per second there began pitching oscillations 
with slowly increasing amplitudes which continued to a speed of about 
620 feet per second for model D and 500 feet per second for model E. 

Over this interval the primary frequencies of oscillation increased from 
U to 7 cycles per second for model E, from 5 to 8 cycles per second for 
model D, and there was a considerable component of the wing first natural 


6 


NACA m L50I29 


bending frequencies present on the wing bending strain-gage traces.. Then 
there occurred an interruption or damping of the oscillations followed 
by final divergent oscillations with frequencies of about 9.0 cycles 
per second for model D and .8.5 cycles per second for model E which 
resulted in failures of the wing after about 7 or 8 cycles. The final 
oscillations were of a much more regular character, approaching a 
divergent sine wave. The speeds at the beginning of the final divergent 
oscillations, 670 feet per second for model D and 580 feet per second 
for model E, are considered to be the- experimental flutter speeds because 
of the regularity of the subsequent oscillations. It may be of interest 
to note that the periods of the low-frequency oscillations for both 
models, from the beginning of a detectable oscillation until wing failure, 
were about 100 chords per cycle, as indicated by the time histories. 

The primary interest in flutter, research centers on the speed and 
. frequency of oscillation. A secondary interest exists in the deflection 
amplitudes and in the phase relationships of the various motions j the 
latter are usually not experimentally determinable with good accuracy. 
From analysis of the oscillograph records of model E over the interval 
t = 5.1*0 to 5.65 seconds during flutter, the amplitudes of wing defor- 
mation and body motions relative to a wing-tip-bending amplitude of 
1 inch were found to average about 0.0065 radian for wing- tip torsion, 
0.020 radian for pitching, and 0.21 inch for vertical translation. The 
wing-tip-bending and wing-tip-torsion amplitudes were determined through 
the assumption of constant lift distribution from wing tip to tip, while 
the amplitude of vertical translation was obtained by integrating the 
normal acceleration twice. The body-pitching amplitudes were determined 
by using the normal acceleration and an assumed lift-curve slope of 
0.080 per degree. The over-all accuracy of the experimentally determined 
amplitudes is thought to be within plus or minus 50 percent, and thus 
they give at least the proper order of magnitude. 

It can be seen from figure 7 that the phase relationship is as 
follows within the limit ±15° of experimental accuracy: at an instant 
when the wing is bent up a maximum, it is also twisted leading edge up 
a maximum and, as interpreted from the normal acceleration, the missile 
is pitched nose up and translated down a maxi mum amount. 


ANALYTICAL INVESTIGATION 
Bases of Analysis 


Only an outline of the analytical treatment used to obtain the 
results will be given here. It is assumed that the theory of linear 
superposition for small disturbances holds. A simplified configuration 
which has four. degrees of freedom in a symmetric oscillation is treated. 
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There are two rigid-body degrees of freedom: pitching about the center 

of gravity, and vertical translation; and two wing-deformation degrees 
of freedom which are approximated by. the first-uncoupled-bending and 
f irst-uncoupled-torsion modal shapes of an ideal uniform cantilever beam. 
For convenience, the aerodynamic coefficients used for a wing section are 
those of incompressible, two-dimensional potential flow and yield forces 
which are proportional to the displacement and motion of that section 
(strip analysis). Spanwise distributions of air forces due to symmetric 
motions of a rigid body are thus constant in magnitude and phasing from 
wing tip to wing tip, including the body intercept. In view of the Mach 
number of the experimental flutter, well below the critical Mach num- 
ber M cr of the airfoil section, aerodynamic coefficients for incompres- 
sible flow are considered satisfactory. Coefficients for compressible 
flow could be substituted, however. 

The configuration has no other horizontal airfoil surface and no 
other aerodynamic forces are assumed on the body. Each half wing is 
u nif orm and- the two half wings, right and left, of each missile are 
treated as being identical. The effect of gravity on an oscillation 
having the frequency of the observed oscillation is considered negligible 
as are the effects of rocket thrust. 


Equations of Equilibrium 

In order to obtain the theoretical conditions for flutter Lagrange's 
equation is employed, as for example in reference 7, to derive four 
equations of equilibrium. With the use of the specified condition of 
harmonic motion the four equations may be written in matrix form 


a ll a 12 a 13 a lh 


~ x l~ 


Wx^l 

a 21 a 22 a 23 a 2U 


x 2 


s fef X2 

a 31 a 32 a 33 a 3U 


x 3 


fipx3 

a iil a U2 a ii3 a UU 


_ x li_ 


fiTx^ 


( 1 ) 


In this equation the quantity 


a ' ^) 2(1 * lg) 


where 


has the 


properties of a structural damping coefficient. (Structural damping 
force is proportional to amplitude and in phase with velocity.) The 
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quantities x-^, x 2 , x^, x^ represent the nondimensional amplitudes 

of oscillation in wing bending, wing torsion, body pitching, and vertical 
translation, respectively. On the basis of the foregoing assumptions, 
for the configuration treated ' 


a ll ~ k ~ A ch 


a 12 - - A oa ) 

a 13 . - Acp ) 

~ 1. £662a^ 

■ 0.6778(2; - Aah ) 


31 ‘ °-33177 ( - ~ - % - A^) 
Is! + ( ~ s + k)x a _ 


21 


r a ■ 

aoo - ■" — A 


22 


a 23 = 1.2732 


aa 

; 2 

a_ 

K 


a 32 ' = 0.5395 

r 2 

a 33 = -7- - Ad p 

• a 3li = ~ A dh 
a Ul = 0.33177a^ 
V = °- 39793a- 


A da 


12 


r a + (-, s * a ) x q _ 

K K 

a 2 |^ = 1.878Ua 2 -L 


l ap 


a li3 ' A cp 

a liii = “ A ch 


where 


A ch - - 1 - f ♦ if 


A ca = A ca (k, a) -.» * g - (i - a)f ♦ i[i ♦ § . (I _ a)fj 
A ah - Aaj-,1 k, a) . - I - (i ♦ a) Aoh 


A aa ' A aaO', a) = - i - a 2 - (i . a)f * (j - a*) 


2G 


i[( 


1 . a 2 )a: . /i „ a \2a 

\ 

2 / k \ 

ll / k \2 ) k 2 
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•A-cp = A ca (k, s ) 

^dh = A a h(k, s) 

^dp = A aa (k, s ) 

^ap ( a > s ) = A aa — (s — s.) A- ah 
A^aik, a, s) = A aa - (s - a)A ca 


and. 

2 

rp Mass moment of missile in pitch about its center of gravity 
K (Span)npb^ 

1 _ Total mass of missile 
K T (Span)npb^ 

The quantity a is the nondimensional (in wing semichords) distance 
from the wing midchord line to the wing elastic axis, positive when the 
elastic axis is behind the midchord. The quantity s is the nondimen- 
sional (in wing semichords) distance from the wing midchord line to the 
center of gravity of the missile, positive when the center of gravity is 
behind the midchord. The quantities F' and G are the real and imagi- 
nary parts, respectively, of the complex function C = C(k) = F(k) + iG(k) 

developed in reference 8. The reduced frequency parameter k = 

It is to be recognized that equation (1) represents a characteristic- 
value problem of aeroelastic harmonic motion. The a^j's depend on the 
reduced frequency k, while ft contains the frequency co. Combinations of 
co and k which result in a specified value of g are the characteristic 
values (eigenvalues) of the system. Knowledge of k and co leads to 
a flutter speed v. The trivial solution x^ = X 2 = x-^ = x^ = 0 is of 
course not sought. 

The quantities p and T are proportional to any mechanical 
restoring force due to body motions in the body-pitching and vertical- 
translation degrees of freedom, respectively; Such restoring forces 
could exist, for example, if the model were mounted on springs in a 
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wind tunnel. When the wing-body combination is of the freely-flying 
type, however, there is no mechanical restoring force; that is, the 
natural (in vacuo) frequencies of the two rigid-body degrees of freedom 
are zero. Hence the characteristic-value solutions for equation (l) 
are sought at the limit as p and T approach zero. 


DISCUSSION AND COMPARISON OF ANALYTICAL AND EXPERIMENTAL RESULTS 
Rayleigh-Ritz Treatment in Four Degrees of Freedom 


.When equation (1) for harmonic motion is treated in the Rayleighs 
Ritz manner there exist in general four roots of flutter speed and fre- 
quency which satisfy the equations of equilibrium, although one or more 
roots may have no physical significance. For the configuration tested 
and analyzed the critical (lowest) flutter speed root (designated root A) 
corresponds to an oscillation involving appreciable proportions of the 
rigid-body stability modes at a frequency which is a fraction of the wing 
natural first bending frequency. The next higher speed root (designated 
root B) corresponds . to the conventional wing flexure- torsion flutter and 
only minute amplitudes of rigid-body motions are present. 

During the flight time after launching several flutter parameters 
were continually changing. Such varying parameters were: air density, 

weight and mass moment in pitch of the missile, and the location of the 
missile center of gravity. These combined changes have only a moderate 
effect, an increase of about 7 or 8 percent, on the calculated flutter 
speed (root A) from t = 0 to the time of wing failure. The calculated 
flutter speeds as functions of time are included in figure 1;. (The 
corresponding changes in flutter frequencies, amplitude ratios, and 
phasings are also quite small.) The intersections of calculated flutter 
speeds and flight speeds give the predicted flutter speeds, 520 feet 
per second for model D and U35 feet per second for model E, as shown in 
figure U. The corresponding predicted flutter frequencies are 6.0 cycles 
per second for model D and U.3 cycles per second for model E. Each con- 
figuration did not flutter, however, until its flight speed was 30 percent 
higher, 670 feet per second for model D and 580 feet per second for 
model E, as also shown in figure U. Thus the analysis was conservative 
in its predictions of flutter speeds. 

With regard to the degree of conservatism of the analytically pre- 
dicted flutter speeds it should be pointed out that even though the 
aspect ratio is 7 (including body intercept), the finite- span correction 
to the flutter speed is quite large for small values of the reduced fre- 
quency parameter k. For models D and E the analytical value of k was 
approximately 0.027 which is quite small for flutter work. Reference 9 
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gives some results of flutter tests made in part to determine the effect 
of finite aspect ratio on the flutter of cantilever wings. The flutter 
was the bending-torsion type and occurred mostly in a range of k from 
about 0.15 to 0.30. According to the results of reference 9 for an 
aspect ratio of 6 to 8, the finite span causes an increase of roughly 
5 to 10 percent in the flutter speed. For such a low value of k as 
was obtained with models D and E it is known that the finite span causes 
.a considerably larger increase in the flutter speed. Therefore, it 
appears that if the proper finite-span corrections were made, an analyti- 
cally predicted flutter speed would agree much better with the experi- 
mentally determined one. Compressibility effects also influence the 
results and can be taken into account. It is of interest to note that 
an oscillation with k = 0.027 has a period of 117 chords per cycle 
since n/k equals the period. This value of the period is well within 
the range of period treated in dynamic- stability work with rigid or 
undeformable aircraft. 

The amplitude ratios and phasing of the various degrees of freedom 
associated with root A for model E when t = 5«2 seconds based on a 
wing- tip-bending amplitude of 1 inch, are as follows: 


Degree of freedom 

Relative amplitude 

Relative phase angle 
(deg) 

Wing-tip bending 

1.0000 inch 

0 

Wing-tip torsion 

0.001*36 radian 

179.1 

Rigid-body pitching 

0.01679 radian 

186.1 

Rigid-body vertical 
translation 

0.21*37 inch 

186.0 


The parameters used were: mass moment of missile in pitch, 13.25 slug- 

feet^; weight of missile, 211 pounds; air density, 0.002260 slug 
per cubic foot, and the parameters of the right wing of model E. The 
wing-tip-torsion amplitude is one-quarter of the body-pitching amplitude, 
and the vertical translation amplitude is one-quarter of the wing-tip- 
bending amplitude. Each degree of freedom is very nearly in phase or 
almost diametrically out of phase with other degrees of freedom, (on the 
basis that positive translation and wing bending are downward, while 
positive angular displacement is leading edge or nose up). Also of 
interest is the virtual identity of phasing of the two rigid-body modes, 
which fact, together with the appropriate amplitude ratios, indicates 
that the missile is pitching effectively about a lateral axis li*.9 inches 
forward of the missile center. of gravity. All the analytical amplitude 
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ratios and phase relationships agree with their experimental counterparts, 
reported under Experimental Results, within the experimental accuracy. 

The amplitude ratios of wing flexure-torsion flutter (root B), based 
on a wing-tip-bending amplitude of 1 inch, are quite minute for the two 
body degrees of freedom, while the wing-tip-torsion amplitude is 
O.O 676 radian with a phasing of about -18°. It should be remarked that 
the ratio of body weight to wing weight ranged from 38 at time 'of 
launching to 31 at time of wing failure for model E; if this ratio were 
smaller, ( more relative body motion might be associated with root B. 

Effect of Various Binary and Ternary Combinations of Freedoms 

Analyses employing various binary and ternary combinations of degrees 
of freedom were made and the results compared with those of the four- 
degree-of -freedom analysis. R'esults are given in table II. The parameters 
used throughout are those of the right wing of model D and those listed at 
the top of table. II. The four-degree-of -freedom analysis is designated 
type (a). . ■" 

Consider root A. The addition of the wing-torsion degree of freedom 
to an analysis not already including it corresponds to a reduction of 
torsional stiffness from an infinite value to the finite experimentally 
determined value. Even though the amplitude of wing-tip torsion in a 
four-degree-of -freedom analysis, type (a), is small compared to body- 
pitching amplitude, addition of wing torsion to analysis type (c) or (e) 
so that they become type (d) or (a.), respectively, causes a significant 
reduction of about 13 percent in the analytically determined flutter 
speed. These reductions indicate that it is potentially dangerous to 
exclude the wing torsion from an analysis of body-freedom flutter. The 
addition of the vertical-translation degree of freedom to analysis 
type (c) or (d), so that analysis type (e) or (a), respectively, is 
obtained, effects an increase of £0 percent in the flutter speed, demon- 
strating clearly that it is not realistic to exclude vertical translation 
from an analysis treating body freedoms, at least for the present con- 
figuration. Thus it seems clear that at least the four degrees of freedom 
employed in this report should be included in an analysis of symmetrical 
body-freedom flutter. The wing-bending stiffness is of predominant 
importance to the flutter speed, root A, since that flutter speed is 
virtually directly proportional to the natural bending frequency, at least 
in the range of low ratio of bending frequency to torsion frequency. 

The flutter speed and frequency of root B, primarily wing bending- 
torsion flutter, are unchanged regardless of which type of analysis, (a), 
(b), or (d), is employed, and for root B the flutter speed is predominantly 
affected by the torsional stiffness in constrast to the predominant effect 
of bending stiffness on the flutter speed of root A. 
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Effect of Relative Positions of Wing and Missile 
Center of Gravity 

Another variable whose effect was studied was the distance from the 
missile center of gravity to the wing leading edge. The flutter parameters 
used were again those of the right wing of model D and those listed in 
table II, other than the fixed distance between the wing leading edge and 
missile center of gravity. The analysis was type (a); that is, it included 
all four' degrees of freedom. 

Figure 8 shows the two calculated roots A and B of flutter speed 
coefficient as a function of the distance of the wing leading edge behind 
the center of gravity of the missile. The experimental flutter speed 
coefficient for model D is included on the figure at a wing position of 
3.38 semichords, even though the experimental parameters differed slightly 
or moderately from those at the top of table II on which the analytical 
curves are based. If the parameters of the analysis coincided with their 
experimental counterparts, the calculated curves would be affected by only 
a few percent. The important point is that the experimental flutter speed 
is about 30 percent higher than the calculated one. It is pertinent to 
recall that the missile has no horizontal airfoil surface other than the 
wing and that the missile would become statically unstable with the wing 
very far forward near the center of gravity. Such a rigid-body static 
instability is not indicated on the figure. Two low-frequency failures 
of -missiles which did have horizontal tails and which had the missile 
center of gravity at. approximately the wing one-eighth chord are reported 
in reference 5 and will be discussed later in this section. 

If analysis type (c) or (d) (see table II) which do not include the 
freedom of vertical translation were employed, it would be predicted, on 
the basis of the aerodynamic forces of two-dimensional incompressible 
potential flow developed in reference 8, that single-degree-of -freedom 
rigid-body pitching instability would be encountered over the approximate 
range of wing position shown in figure 8. This single-degree pitching 
instability was investigated analytically by Smilg in reference 10. The 
quantity I<x/npb^ of reference 10 (where I a represents the mass moment 
in pitch of the missile divided by the span) is about ten to eleven 
thousand for models D and E. The predicted borderline flutter speed 
(root A), if no vertical translation were permitted, would be zero in the 
critical range of wing location, since no mechanical restraint of pitching 
exists; that is, the in vacuo natural rigid-body pitching frequency is 
zero. It appears, however, that such an unhappy possibility is precluded, 
even in two-dimensional incompressible flow, by the actual existing free- 
dom of vertical translation of the missile in flight, at least for the 
configuration treated. 
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The Introduction includes mention of low-frequency failures of two 
models, 6 and 7, reported in reference 5. Each of these models had a 
pair of unswept rectangular uniform model wings 2k by 8 inches. The 
wings of model 6 were designed to be similar to, though slightly smaller 
than, the wings of models I) and E of the present report, and the wing- 
first-bending and first-torsion frequencies of model 6 were practically 
identical to those of model D. The leading edges of the wings of models 6 
and 7 were located slightly ahead of the missile center of gravity, 

-0.25 semichord for model -6, and -0.38 semichord for model 7. The com- 
plete models 6 and 7 had weights and mass moments in pitch approximately 
one-third as great as those of models D and E. Each model included a 
horizontal tail having an area equal to about 20 percent of the wing 
area, and a tail length of about one-quarter of the total tip-to-tip 
wing span. 

Thus it may be seen that models 6 and 7, except for the added hori- 
zontal tails, were dynamically fairly similar to a hypothetical model D 
or E which has had its wings moved forward near to the missile center of 
gravity. In the light of such observation, it is interesting to examine 
figure 8 for a small negative value of the abscissa or wing position. 

In that region the speed of the flutter involving body motions is pre- 
dicted to be lower than the speed for wing bending-torsion flutter of a 
modified model D or E. 

All these circumstances lead to the conjecture that the horizontal- 
tail volumes of models 6 and 7 were sufficiently .large to prevent a 
static instability, but not sufficient, or the tails were improperly 
placed, to prevent a low-frequency body-freedom flutter. 


CONCLUSIONS 


The following conclusions may be drawn from flight tests and a 
mathematical analysis made to investigate a type of flutter involving 
wing deformation and body motions of a particular configuration and 
certain modifications of that configuration: 

1. The possibility of / a flutter type of dynamic instability 
involving appreciable proportions . of rigid-body motions as well as wing 
structural deformations has been predicted analytically and confirmed 
experimentally. The flutter was entirely different from conventional 
wing bending-torsion flutter in that the frequency was about one-half 
the wing natural first-bending frequency. The period was approximately 
100 chords per cycle which is well within the range of period found in 
rigid-body control-free stability work. 
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2. A flutter analysis based on two-dimensional incompressible flow 
provided a conservative (by 30 percent) prediction of the airspeed at 
which the flutter would occur for the configuration studied. 

3. Wing-bending stiffness rather than wing-torsional stiffness has 
the predominant effect on speed of the body-freedom flutter. 

U. At least four degrees of freedom, wing bending, wing torsion, 
body pitching, and body vertical translation, should be included in an 
analysis of symmetrical flutter involving wing deformations and body 
motions. 

5. Analysis predicted that if the wing were moved rearward from a 
position of coincidence of wing quarter chord and missile center of 
gravity (no horizontal tail present) the speed of flutter involving body 
motions would first decrease sharply, until a certain wing position 
relative to the center-of -gravity location was reached, then increase 
slowly, as the wing was moved farther to the rear. 


Langley Aeronautical Laboratory- 

National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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TABLE I 


WING PARAMETERS 


Parameters 

Model D 

Model E 

Left wing 

Right wing 

Left wing 

Right wing 

NACA section . . . 

65A006 

65A006 

65A006 

65a 006 

^cr 

0.81* 

0.81* 

0.81* 

0.81* 

1 

29.875 

29.875 

30 

30 

Act 

3 

3 

3 

3 

b 

0.1*19 

0.1*19 

0.1*17 

0.1*17 

x ce • 

38.8 

37.6 

1*0.0 

1*0.91* 

V D 

Xea 

33.5 

35.2 

36.75 

37.5 

a 

-0.33 

-0.296 

-0.265 

-0.25 

a + Xa , 

-0.221* 

-0.21*8 

-0.20 

-0.1812 

1 

21*. 6 

25.1 

31.11 

31.08 

k (stnd.) 





r 2 

ra 

0.2255 

■ 0.2105 

0.183 

0.1837 

% 

19.5 

20 

16 

16.5 

fh? 

115.0 

115.5 

99 

106 

ft 

131* 

138.5 

118.5 

119 

fa 

131 

137.5 

117.0 

117.5 

GJ . 

3.36 > 10^ 

3.61 x io£ 

2.552 x 10? 

2.555 x 10^ 

El 

2.76 x io^ 

2.93 x 10$ 

2.31 x 10^ 

2.31* x 10^ 

gh 

0.05 


o.ol* 

o.ol* 

6a 

0.01 

0.02 


0.01 


'\naca3^’ 
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TABLE II 


MISSILE PARAMETERS 


[Missile center of gravity is 17 inches ahead of wing 
leading edge: missile mass moment in pitch is 
lii.8 slug-ft^j missile weight is 236 lb; air 
density p is 0.002272 slug per cu ft7| 


Type of 
analysis 

Degrees of 
freedom included 

(1) 

Flutter-spee'd 

coefficient 

v/bcoh 

Flutter 

frequency ratio 

co/« h 

Root A 

( 2 ) 

Root B 
(3) 

Root A 

Root B 

a 

F 

T 

p 

V 

9.6 

22.5 

0.26 

3. ho 

b 

F 

D 

- 

B 

— 

22.5 

— 

3. h0 

c 

F 

1 

■ .p 


7.3 

— 

0.20 

— 

d 

F 


p 

B 

6.U 

22.5 

0.17 

3.1*0 

e 

F 

■ 

-D 

j. 

D 

11.2 

— 

0.287 

— 


^Designation of degrees of freedom 
F - wing flexure (first bending) 

T - wing torsion (first torsion) 

P - rigid-body pitching 
V - rigid-body vertical translation 

O 

Root A - critical body-freedom root 
%oot B - critical wing bending-torsion root 
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Figure'!.- Layout of flutter test vehicle, model E. All dimensions are 

in inches . 
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Figure 2.- Variations of model E parameters with time. 










Figure 3 *- Model E in launching position. 
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Figure 4.- Velocity plotted against time. 
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Figure 5*- Air density and Mach number variations with time. 
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Figure 6.- Variations of longitudinal acceleration with time. 
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Figure 7.- Portion of telemeter record of model E. 
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